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Abstract 

We apply a framework for the description of random tilings without height representation, 
which was proposed recently, to the special case of quasicrystalline random tilings. Several 
important examples are discussed, thereby demonstrating the consistency of this alternative 
description with the conventional one. We also clarify the latter by deriving a group theoretic 
criterion for the validity of the first random tiling hypothesis. 

1 Introduction 

The systematic study of random tilings arose from the insight that they may serve as models 
for entropically stabilized quasicrystals [^|, 0|. As such, they are an important and interesting 
alternative to the perfect tilings based upon the projection method. Although random tilings were 
used intensively to study transport properties and dynamical properties of real quasicrystals, not 
much effort was spent to deduce random tiling properties from general grounds || |^] . We mention 
the two random tiling hypotheses stated by Henley in 1991, which serve as a starting point to 
infer diffraction properties, using ideas of the Landau theory of phase transitions. Recently, a more 
general approach was investigated |l6L [H], ^ , mainly for two reasons: On the one hand, the concept 
of entropic stabilization applies also to crystalline solids which therefore should be included in a 
more general description. On the other hand, the quasicrystalline random tilings share a special 
feature which stems from the fact that these tilings are derived from the perfect ones: They have 
a height representation - each tiling can be embedded as a surface in a higherdimensional space. 
Since tiling surfaces with equal (mean) density of prototiles have equal (mean) slope, the slope 
parameters can be used to describe the random tiling ensemble. This is the so-called phason strain 
parametrization which governs the conventional description. 

This description has a number of shortcomings: First of all, it only allows the description of 
those random tilings which do allow a height representation - being the vanishing (though interest- 
ing) minority of all possible tilings. Second, symmetry analysis in this framework is constrained to 
geometric symmetries, whereas there may be other relevant symmetries of the tiling ensemble. An 
example are colour symmetries modelling chemical disorder, which is shown to play an important 
role for real quasicrystals |^|, [l(J. Third, it was expected that a more general analysis could also 
clarify the origin of the random tiling hypotheses mentioned above. 

This was the viewpoint which led to the random tiling description proposed in ]l^] . There, the 
grand-canonical tiling ensemble was considered where prototiles are energetically degenerate, and 
the chemical potentials of the different prototiles or their (mean) densities are the only macroscopic 
parameters. It was explained how to perform a symmetry analysis in this framework. This led 
to a proof of a generalization of the first random tiling hypothesis which asserts that the point 
of maximum entropy is a point of maximum symmetry. In addition, the range of validity of 
the second hypothesis could be analyzed using the grand-canonical setup, and an exactly solved 
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(crystallographic) counterexample could be given. A rigorous treatment of diffraction is also 
possible using this setup 

The examples presented so far were taken only from the crystallographic case. It remained 
to apply the proposed description scheme to quasicrystallographic random tilings and to show 
the consistency of the two approaches. This is where this paper aims at: The next two chapters 
are devoted to the description of two important examples of quasicrystallogrphic random tilings, 
the Ammann-Beenker tiling |Q and the square-triangle tiling Jl5|, [l^, [ll], |l2). We perform a 
symmetry analysis along the lines of ]l6| and compare the approach with the conventional one, 
thereby showing the consistency of the two approaches. In the Appendix, we shortly review and 
clarify the conventional description. We analyze the connection between maximal symmetry and 
vanishing phason strain and derive a group theoretic criterion. If the criterion is valid - which 
is the case in all examples we met - the first random tiling hypothesis is fulfilled in its original 
formulation: The point of maximum entropy is a point of maximum symmetry [fl6| , and maximum 
symmetry in turn implies vanishing phason strain. A comparing discussion of the two approaches 
concludes the paper. 

2 The Ammann-Beenker Random Tiling 

The prototiles of the Ammann-Beenker random tiling are two squares and four 45°-rhombi as 
shown in Fig. |[ There is no further matching rule apart from the face-to-face tiling condition, 
which imposes the relation pj + y] = 1 on the prototile densities. Tilings can alternatively 




Figure 1: Prototiles of the Ammann-Beenker random tiling 



be viewed as world lines of particles of two different kinds, as is indicated by the decoration. 
It was shown in |I3] that there is a nonlinear constraint among the different prototile densities 
which reduces the number of independent variables to four. This constraint was obtained using 
the height representation of this model, which is given below. It should be noted that it can be 
independently derived relating the tile description to the line interpretation. This method was 
introduced in ||. The constraint reads 

pip 2 = 2 (ctict 3 + 0-20-4) ■ (1) 

The derivation imposes periodic boundary conditions. It is highly plausible that the constraint 
holds (asymptotically) also in the case of free boundary conditions, though there is no proof of this 
assumption. Complementary, the line representation suggests four independent parameters since 
it is possible to fix the densities of different lines as well as their mean direction independently. 
Moreover, care has to be taken in choosing a set of independent parameters: The four rhombi 
densities, for example, allow the determination of the square densities only up to a permutation, 
due to the nonlinear constraint. This subtlety does not arise if, for example, three rhombi densities 
and one square density are chosen. The quadratic invariants given below are insensitive to this 
because only absolute values of square density differences occur. 

At the point of maximum symmetry, the four rhombi densities and the two square densities 
have the same value. The nonlinear constraint fixes them to 

Pi = 1/4, Ui = 1/8. (2) 
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This is the point where the squares and the rhombi each occupy half of the tiled area. We 
determine the second order expansion of the entropy density according to 2?8-symmetry. We 
introduce reduced prototile densities 

n = Pl - 1/4, Sl = tr, - 1/8, (3) 

such that the point of maximum symmetry is the origin. The vector space relevant for the 

symmetry analysis is obtained through linearization of the constraints. It is four-dimensional 
and determined by 

n + r 2 = 0, si + s 2 + s 3 + s 4 = 0. (4) 

The analysis of the symmetries in the vector space of the independent, reduced prototile densities 
leads to a second order entropy expansion of the form 

s (r,s) = So -i^A i /W(r )S ) + ... (5) 

i 

There are three (positive) elastic constants with invariants 

J«(r,«) = 4r?, (6) 
I^(r,s) = (s 1 + s 3 ) 2 , 

I i3) (r )S ) = \(s!-s 3 ) 2 + ( Sl+ 2 S2 + s 3 ) 2 . 

This can be obtained by decomposing the representation of the symmetry group in the vector 
space of the reduced, independent density parameters according to |]l6|| , in full analogy to the 
height representation treatment which is presented in the Appendix. 

We now derive the relation to the elastic constants defined from the height representation 
approach. This is done by expressing the expansion (^) in terms of the slope parameters of the 
height representation. The positions of tiling vertices are integer linear combinations of vectors of 
the regular 8-star. They therefore belong, viewed as complex numbers, to the module Z[£] with 
£ = e 2 ™/ 8 . We characterize the height function algebraically. To this end, we recall some facts 
concerning the associated cyclotomic field [jl3|. £ is a root of the eighth cyclotomic polynomial 
Pg(x) = x 4 + 1. The other roots are £ 3 , £ 5 = — £ and £ 7 = £. The automorphism group on this set 
of primitive roots of unity is the Galois group Q ~ C 2 x C 2 . The automorphisms can be uniquely 
extended to automorphisms of the corresponding module. The usual height function h on the set 
of possible vertex positions is defined by 

h(a) = a{a) (7) 

with a E Q 3 a : £ i— > £ 5 . The set of all embedded vertex positions constitutes the four- 

dimensional primitive hypercubic lattice Z 4 , and each tiling corresponds to a two-dimensional 
surface. This is the geometric origin of the height representation. It is possible to define a height 
function on the prototiles by linear extrapolation of the height function on the corresponding 
vertices. The new macroscopic parameters to describe the tiling ensemble are the components of 
the phason strain tensor 



E 2 ,l -©2,2 

in cartesian coordinates. The phason strain is defined as (mean) slope of the height function. This 
leads to the relations 

E\.i = pi -pi + (oi + 02) - (03 +04), (8) 
E 2 . 2 = pi - p 2 - (<Tl + o 2 ) + ((T 3 + 04), 
Ei f2 = 2(0-4 - cr 3 ), 
Eh,i = 2(0-1-02). 
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We observe E = at the point of maximal symmetry, which already follows from group theoretical 
considerations given in the Appendix. In order to write the entropy expansion (|J) in terms of 
the phason strain, the relation between phason strain and independent tile densities has to be 
invertible at the point of maximum symmetry. This is not the case if the four rhombi are taken 
as parameters. If three rhombi densities and one square densities are taken instead, the relation 
(||) is invertible in the whole phase space. This leads to an expansion 

a(E) = s ^^K 1 (E hl +E 2 , 2 ) 2 ~^K 2 (E h2 -E 2A ) 2 (9) 

~K 3 ((E hl - E 2 , 2 f + (E h2 + E 2A ) 2 ) + ... 

This expression is known from the literature ||l4| and is rederived in the Appendix using the 
symmetry analysis in the height representation. The relation between the elastic constants is 

Ai = 4RTi > 0, A 2 = 16K 2 > 0, A 3 = 16K 3 > 0. (10) 



3 The Square- Triangle Random Tiling 

The prototiles of the square-triangle random tiling are three squares and four equilateral triangles 
as shown in Fig. ^|. There is no further matching rule apart from the face-to-face tiling condition, 
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Figure 2: Prototiles of the square-triangle random tiling 



which imposes the relation Pi a i = 1 on the prototile densities. Tilings can alternatively be 
viewed as world lines of particles of two different kinds, as is indicated by the decoration. Using 
this representation, it is easy to see that triangles always occur in pairs. Between the (mean) 
prototile densities we therefore have the relations p x = p 3 and p 2 = p 4 . It was shown in ||] that 
there is a nonlinear constraint among the different prototile densities which reduces the number 
of independent variables to three. The constraint reads 

16pi/92 = 3 (CTi<72 + 0-203 + O3CT1) . (11) 

It should be mentioned that the derivation imposes periodic boundary conditions. The same 
constraint is believed to hold (asymptotically) in the case of free boundary conditions. This is 
suggested from the complementary observation that the line representation allows three indepen- 
dent parameters: the densities of the lines of different type and the frequency of one type of 
crossing. Due to the nonlinear constraint, care has to be taken in choosing a set of independent 
parameters: The three square densities, for example, allow the determination of the triangle densi- 
ties only up to a permutation. This subtlety does not arise if two square densities and one triangle 
density are chosen instead. The quadratic invariants given below are insensitive to this because 
only absolute values of triangle density differences occur. 

At the point of maximum symmetry, the four triangle densities and the three square densities 
have the same value. The nonlinear constraint fixes their values to 



pi = 1/8, ^ = 1/6. (12) 

This corresponds to the situation where the squares (resp. triangles) occupy half of the tiled 
area. We determine the second order expansion of the entropy according to 2?i2-symmetry. We 
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introduce reduced prototile densities 



n = pi -1/8, Si = (Ti- 1/6, (13) 

such that the point of maximum symmetry is the origin. The vector space relevant for the 
symmetry analysis is obtained through linearization of the constraints. It is three-dimensional 
and determined by 

n + r 2 = 0, si + s 2 + s 3 = 0. (14) 

The analysis of the symmetries in the vector space of the independent, reduced prototile densities 
leads to an entropy expansion of the following form, up to second order: 

s(r,s) = ao - i£ Ai/W(r, «) + ... (15) 

i 

There are two (positive) elastic constants with invariants 

7 (r) (r,s) = Ar\, (16) 
&Hr,s) = 6(sl- Sl S2 + si). 

The invariants are obtained from the decomposition of the representation of the symmetry group 
into its irreducible components, according to We now derive the relation to the elastic 

constants of the height representation approach. The positions of tiling vertices are integer linear 
combinations of vectors of the regular 12-star. They therefore belong, viewed as complex numbers, 
to the module Z[£] with £ = e 27T ^ 12 . We characterize the height function algebraically. Therefore, 
we remind of some facts concerning the associated cyclotomic field pjfl . £ is a primitive root of 
unity, annihilating the 12th cyclotomic polynomial Pi 2 (x) = x 4 — ar -f 1. The other roots are 
£ 5 , £ 7 = — £ and £ n = £. The automorphism group on this set of primitive roots of unity is the 
Galois group Q ~ C 2 x C2. The automorphisms can be uniquely extended to automorphisms of 
the corresponding module. The usual height function h on the set of possible vertex positions is 
defined by 

h(a) = a(a) (17) 

with a 6 Q 9 a : £ 1— > £ 7 . It is possible to define a height function on the prototiles by 

linear extrapolation of the height function on the corresponding vertices. The new macroscopic 
parameters to describe the tiling ensemble are the components of the phason strain tensor 

E — I ^1,1 -Ei, 2 

-E-2,1 E2.2 

in cartesian coordinates. The phason strain is defined as (mean) slope of the height function. This 
leads to the relations 

Ex,x = ax - i(cr 2 + 0-3) + (pi + Pz) - (P2 + Pa), (18) 

#2,2 = -fl + -(^2 + ^3) + {Pi + A3) - (p2 + Pi), 

^l, 2 = ^2,1 = -— (0-3 — a%). 

We find E = at the point of maximal symmetry, which already follows from group theoretical 
considerations given in the Appendix. Since the model has three degrees of freedom, as follows 
from the density analysis, we have to expect a special property of the height function which reduces 
the number of independent slope parameters, in contrast to the previous example. It is indeed 
easily seen that each height function has the irrotationality property fig] 



3? (^<j)h{z)dz^ = 0. 



(19) 
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The equality of the off-diagonal elements is a consequence of the irrotationality property (|19|) . (It 
is also possible to define a height function by Q 3 a : £ i— > £ 5 . In this case, the irrotationality 
property reads 5ft (§ h(z) dz) = 0.) 

In order to write the entropy expansion ( |l5| ) in terms of the phason strain, the relation between 
phason strain and independent tile densities has to be invertible, at least at the point of maximum 
symmetry. If the three squares are taken as parameters, this is not the case in the phase with 
C4-symmetry. If two square densities and one triangle densities are taken instead, the relation ( |l8| ) 
is invertible in the whole phase space. This leads to an expansion 



which is known from the literature fllq] and is rederived in the Appendix using the symmetry 
analysis in the height representation. The relation between the elastic constants is 



4 Conclusion 

We performed a symmetry analysis of well-known planar quasicrystalline random tilings using the 
prototile density approach. We were able to show the consistency to the conventional approach via 
the height representation of these models. For the square-triangle tiling, the approach illuminated 
the irrotationality property whose appearance was rather unexpected before. Whereas it seems to 
be natural to describe random tilings using the prototile densities as macroscopic parameters, the 
approach is commonly hard to follow because it requires the knowledge of all (possibly nonlinear) 
constraints in order not to overestimate the number of elastic constants. The height representation, 
in turn, has the disadvantage that is is not obvious whether the slope parameters are independent 
and whether they constitute a sufficient set to describe all degrees of freedom of the ensemble. 
Furthermore, to describe other than geometric symmetries, one has to derive the relation between 
phason strain and the prototile densities explicitly. Therefore, the two approaches are in fact 
complementary in deriving properties of random tilings. The exactly solved eightfold random 
tiling, for example, must possess a constraint additional to the ones derived in J5j, since the height 
representation and the density description yield different numbers of elastic constants otherwise, 
as can be shown. Research in this direction is in progress. 

Acknowledgments 

The author is grateful for discussions with Michael Baake, Jan de Gier, Joachim Hcrmisson and 
Moritz Hoffc. He also thanks Michael Baake for several useful comments on the manuscript. 
Financial support from the German Science Foundation (DFG) is gratefully acknowledged. 

Appendix: Symmetry and Phason Strain 

We shortly expose the symmetry analysis of random tilings which permit a height representation, 
in the framework of our viewpoint exposed elsewhere Jl6| | . In particular, this will lead to a criterion 
for the validity of the first random tiling hypothesis [M which states that the point of maximum 
entropy occurs at vanishing phason strain. 

Let a random tiling ensemble of the d|| -dimensional vector space V|| be given. A height rep- 
resentation assigns to each tiling a d|| -dimensional surface in a higher-dimensional vector space 
V = Vji ® V± with the property that surfaces of tilings with equal prototile (mean) density have 
the same (mean) slope in V. In this situation, the (mean) surface slope can be used to parametrize 
the grand-canonical tiling ensemble. We discuss the influence of geometric symmetries on this de- 
scription. The (mean) tiling slope can be represented by a linear map E : V\\ —> V± , usually called 



s(E) = s - -K» (tr(£)) 2 + -K e det(E) + . . . , 



(20) 




(21) 
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phason strain. Symmetries are bijections on the (mean) prototile densities which leave the entropy 
density invariant. Let Du denote the representation of the geometric symmetries in the physical 
space Vji . The height representation induces a representation D± of the geometric symmetries in 
the internal space Vj_ . The phason strain tensor transforms under geometric symmetries according 
to 

E = D^EDT, 1 . (22) 

We now focus on the point of maximum symmetry. E is invariant at this point, since this is by 
definition true of the (mean) prototile densities, which determine E uniquely. If the representations 
D|| and Dj_ are irreducible and not equivalent, Schur's lemma yields E = 0. In the general case 
we have 

• The phason strain tensor vanishes at maximal symmetry if no irrep in internal space is 
equivalent to an irrep in physical space. 

If this criterion is fulfilled, the first random tiling hypothesis in the formulation of Henley is 
satisfied since the point of maximum entropy is always a point of maximum symmetry, as was 
shown in ]jq| . On the other hand, the criterion is fulfilled in all situations we met due to the 
special construction of the height function. 

Symmetry analysis may alternatively be performed in the tensor product Vy <8> Vj_, which is 
most advantageous when discussing the influence of symmetries on the entropy density. In this 
case, E is regarded as a vector with components of the matrix E : Vji — * V±. The geometric 
symmetries are represented via 

E = D-E= (jo^ 1 ) ' ®-°-lJ (23) 

Since the point of maximum symmetry is by definition a fixed point, we conclude: 

• The phason strain at the point of maximum symmetry consists of components in direction 
to the trivial ld-irreps of the representation in Vj| ® V±. 

• In particular, the phason strain vanishes at maximal symmetry, if the representation in 
Vji <8> Vj_ does not contain trivial parts. 

This criterion is equivalent to the statement given above, as follows from a closer look at the 
Clebsch-Gordan decomposition of the product representation. 

We give an expansion of the entropy density of the random tiling ensemble in terms of the 
phason strain at E — 0. Without loss of generality we assume that the symmetries are represented 
by orthogonal mappings in Vji ® V±. The entropy density is invariant under symmetry transfor- 
mations, s(E) = s(E). Since E = is a fixed point of the action of the symmetry group, the 
symmetries commute with the Hessian H , the second derivative of the entropy density, 

H\ E=0 D = DH\ E=Q . (24) 

According to Schur's lemma, H\ E=0 acts trivially on the irreducible components in Vji (8> Vj_. 
Therefore, H\ E _ is a linear combination of projectors onto the irreducible components, 

tf| £ =o = -E A * p(l) - ( 25 ) 

i 

The (positive) Xi are the elastic constants of the random tiling, arising from the height represen- 
tation. The second order term of the entropy expansion is of the form 

s 2 (E,E) = -lj2 X i l{i) ( E > E )> (26) 
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where jW {E, E) denotes the quadratic invariants 

I®(E,E) = (E,P®E) = \\P^E\\ 2 . (27) 

(-, •) denotes the scalar product in Vji ® V± and || • || the induced norm. 

We treat as example two-dimensional random tilings with dihedral symmetry T> n . This com- 
prises all models discussed above. In the physical space, generators of the symmetries may be 
given by a rotation Ru and a reflection S\\ . We assume that the height representation induces 
a representation in the internal space with generators R± = —R\\ and S± — S\\ . (This is the 
canonical choice for n even.) Since the corresponding representations Du and D± are irreducible 
and not equivalent, we immediately conclude E = at maximal symmetry. Let us determine the 
quadratic invariants. The representation of the symmetry group in Vji ® V± decomposes into three 
parts, a two-dimensional rotation and two one-dimensional reflections. If we introduce standard 
bases, this yields three elastic constants with invariants 

lM(E,E) = i(£ M +£ 2 , 2 ) 2 , (28) 

I^\E,E) = i(£ x , 2 - £ 2>1 ) 2 , 

I^\E,E) = ^(E hl -E 2 , 2 ) 2 + ^(E 1 . 2 + E 2 ,i) 2 . 

Mind that the square-triangle tiling is described by only two elastic constants, since the off-diagonal 
elements in E are equal, due to the irrotationality property (|l9|). 
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